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Abstract 

We have made a precise measurement of the absolute branching fractions of B® — > Di*^ + 
decays using 121.4 fb _1 of data recorded by the Belle experiment running on the T(5S) resonance. 
The results are B{B° S ^D+Dj) = (0.58 ±om ±0.13)%, B{B° S D*±Df) = (1.8 ±0.2 ±0.4)%, 
and B(B° S D* s + D*~) = (2.0 ± 0.3 ± 0.5)%; the sum is B(B° S -> D$* )+ dP~) = (4.3 ± 0.4 ± 
1.0)%. Assuming B^ D ( s * )+ D^' saturates decays to CP -even final states, the branching fraction 
constrains the ratio Ar s / cos <j) s , where Ar s is the difference in widths between the two B s -B s mass 
eigenstates, and <j) s is the CP-violating phase in B s -B s mixing. For <j) s = we obtain Ar s /T,, = 
0.090 ± 0.009 (stat.) ± 0.023 (syst.), where T s is the mean decay width. We also measure for the 
first time the longitudinal polarization fraction of B® ^ D* + D*~ ; the result is 0.12 1qo5- 
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Decays of B s mesons help elucidate the weak Cabibbo-Kobayashi-Maskawa structure of 
the Standard Model (SM). B s decays can be studied at e + e~ colliders by running at the 
T(5S f ) resonance, which subsequently decays to B^B^ pairs. We have used this method 
previously [lj to study B® — > D (*) + Z)(*)~ decays using 23.6 fb" 1 of data. Here we substantially 
improve this measurement using 121.4 fb _1 of data. In addition to the five-times-larger data 
set, there are other improvements to the analysis: the data have been fully reprocessed 
using reconstruction algorithms with higher efficiency for 7r°'s and low momentum tracks; 
we use larger control samples to evaluate systematic uncertainties; and we take background 
probability density functions directly from data rather than from simulation. We also make 
the first measurement of the fraction of longitudinal polarization (f L ) of B®^-D* + D*~. 

As in our previous study, we reconstruct the final states D+D~, D* + D~ + D*~D+ (= 
D^DJ), and D* + D*~. These are expected to be mostly CP-even, and their partial widths 
are expected to dominate the difference in widths between the two B s -B s CP eigenstates, 
Arf p jsj. This parameter equals Ar s / cos0 s , where Ar s is the decay width difference 
between the mass eigenstates, and <p s is the CP- violating phase in B s -B s mixing jsj]. Thus 



the branching fraction gives a constraint in the Ar 
can receive contributions from new physics (NP) 



parameter space. Both parameters 
51. Previous constraints on Ar o and 



NP contributions to <p s were obtained from a time-dependent angular analysis of B s — > J/ip 



± 
9]. 



decays [6 

At the T(5S) resonance, the e + e~ — > bb cross section is measured to be = 0.34C 
0.016 nb, and the fraction of T(5S) decays producing B s mesons is f s = 0.171 ± 0.030 
Thus the total number of B S B S pairs produced is N , = (121.4 ftT 1 ) ■ a bl ■ f s = (7.11 ± 

s s 

1.30) x 10 6 . Three production modes are kinematically allowed: B S B S , B S B* or B*B S , and 
B*B*. In this analysis we obtain the signal yield for the last mode, for which the fraction 



(/d»-d«) is 0.870 ± 0.017 [lOj. The B* decays via B* a — > £? c 7, and the 7 is not reconstructed. 



'B*B* 

The Belle detector running at the KEKB e + e~ collider 11] is described in Ref. 12[. For 
charged hadron identification, a likelihood ratio is formed based on dE/dx measured in the 
central tracker and the response of aerogel threshold Cerenkov counters and time-of-flight 
scintillation counters. A likelihood requirement is used to identify charged kaons. This 
requirement is 86% efficient for K ± and has a 7r ± misidentification rate of 8%. 

We reconstruct B° S ^D+D;, D^D^aad D* s + D* s ~ decays in which D+^4>n + , K° S K + , 



K*°K+, (j)p + , K° S K*+, and K*°K*+ 



131 ] . Neutral Kg candidates are reconstructed from 
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7r + 7T - pairs having an invariant mass within 10 MeV/c 2 of the nominal K s mass [14j and sat- 
isfying momentum-dependent vertex requirements. Charged tracks are required to originate 
from near the e + e~ interaction region and, with the exception of tracks from K s decays, 
have a momentum p > 100 MeV/c. Neutral K*° (charged K* + ) candidates are reconstructed 
from a K + ir~ (K s ir + ) pair having an invariant mass within 50 MeV/c 2 of m K ,. Candidate 
4> mesons are reconstructed from K + K~ pairs having an invariant mass within 12 MeV/c 2 
of m^. Charged p + candidates are reconstructed from 7r + 7r° pairs having an invariant mass 
within 100 MeV/c 2 of m p+ . The tt° candidates are reconstructed from 77 pairs having an 
invariant mass within 15 MeV/c 2 of m n0 , and with each 7 having an energy E > 100 MeV. 

The invariant mass windows used for the reconstructed candidate (denoted Df) are: 
±10 MeV/c 2 (~3cr) for the three final states containing K* candidates, ±20 MeV/c 2 (2.8<r) 
for <pp + , and ±15 MeV/c 2 (~ 4a) for the remaining two modes. For the three vector- 
pseudoscalar final states we require |cos# hel | > 0.20, where # hel is the angle between the 
momentum of the charged daughter of the vector particle and the direction opposite the 
momentum, evaluated in the rest frame of the vector particle. 

We combine candidates with photon candidates to reconstruct D* + — > -D+7 decays. 
We require E > 50 MeV in the e + e~ center-of-mass (CM) system, and that the energy 
deposited in the central 3x3 array of cells of the electromagnetic cluster exceeds 85% of 
that deposited in the central 5x5 array of cells. The mass difference M~ + — * s 
required to be within 12.0 MeV/c 2 of the nominal value. This requirement (and the 
mass windows) are determined by optimizing a figure-of-merit S/y/S + B, where S is the 
expected signal based on Monte Carlo (MC) simulation and B is the background estimated 
from either MC or mass sideband data. 

We select B® — > D+D~, D^DJ, and D* + D*~ decays using two quantities evaluated 
in the CM frame: the beam-energy-constrained mass M bc = \J E^ jC&m — p 2 B , and the energy 
difference AE = E B — E hcam , where p B and E B are the reconstructed momentum and energy 
of the Bg candidate, and -Ebeam is the beam energy. We determine signal yields by fitting 
events satisfying 5.25 GeV/c 2 < M bc < 5.45 GeV/c 2 and -0.15 GeV < AE < 0.10 GeV. 
Because the 7 from B* — > B s -f is not reconstructed, the modes T(55) — >• B S B S , B S B* and 
B*B* are well-separated in M bc and AE. We expect only small contributions from B S B S 
and B S B* events and fix these contributions relative to B*B* according to our measurement 
using B® — > D~ix + decays We quote fitted signal yields from B*B* only and use these 
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to determine the branching fractions. 

Approximately half of selected events have multiple B® — > candidates. These 

typically arise from photons produced via ir° — > 77 that are wrongly assigned as D* daughters. 
For these events we select the candidate that minimizes the quantity 



1 



(2 + N) 



E 



M M n 

D s D > 



-, 2 



E 

D* 



AM - AM 



-,2, 



a 



AM 



where AM = M~, — M~, and AM = M_,+ — M^ + . The summations run over the two 
daughters and the N (=0,1,2) D* + daughters of a B° s candidate. The mean masses 
M ( ») and widths a M and <r AM are obtained from MC simulation and calibrated for data- 
MC differences using a large B° — >■ D~ control sample from T(45 l ) data. According to 
MC, this criterion selects the correct candidate 83%, 73%, and 69% of the time for D + D~, 
D^DJ, and D* + D*~ states, respectively. 

We reject background from e + e~ — >qq (q = u, d, s, c) events using a Fisher discriminant 



based on a set of modified Fox- Wolfram moments [15|. This discriminant distinguishes jet- 
like qq events from more spherical B^B^ events. With this discriminant we calculate 
likelihoods L s and C q g for an event assuming the event is signal or qq background; we then 
require C S /(C S + C q q) > 0.20. This selection is 93% efficient for signal and removes more 
than 62% of qq background. 

The remaining background consists of T(5S)^BWbW->D+X, T(5S)^BBX (i.e., bb 
hadronizes to B°, B°, or 5 ± ), and B s -> Dfj(23 17) £>W, £>^(2460)L>W, or DfDjir . The 
last three processes peak at negative AE, and their yields are estimated to be small using 
analogous B d —^D s jD^ branching fractions. We thus consider them only when evaluating 
systematic uncertainty due to backgrounds. All selection criteria are finalized before looking 
at events in the signal region. 

We measure signal yields by performing a two-dimensional unbinned maximum-likelihood 
fit to the M hc -AE distributions. For each sample, we include probability density functions 
(PDFs) for signal and qq, BWbW -»■ D+X, and T(5S) -+ BBX backgrounds. As the 
backgrounds have similar M bc , AE shapes, we use a single PDF for them, taken to be 



an ARGUS function [16j] for M hc and a first-order Chebyshev function for AE. The two 
parameters of the Chebyshev function are taken from data in which one of the D + candidates 
is required to be within the mass sideband. 

The signal PDFs have three components: correctly reconstructed (CR) decays; "wrong 



combination" (WC) decays in which a non-signal track or 7 is included in place of a true 
daughter track or 7; and "cross- feed" (CF) decays in which a D^DJ (D* + D*~) is recon- 
structed as a D+D~ (D+D~ or D^DJ), or a D+D; (D^DJ) is reconstructed as a D^DJ 
or D* + D*~ (D* + D*~). In the former case, the 7 from £>*+— )-.D+7 is lost and AE is shifted 
down by 100 — 150 MeV; this is called "CF-down." In the latter case, an extraneous 7 is 
included and AE is shifted up by a similar amount; this is called "CF-up." In both cases 
M bc remains almost unchanged. 

All signal shape parameters are taken from MC and calibrated using B° s — >D^~ir + and 
B° — > D~ decays. The CR PDF is taken to be a Gaussian for M bc and a double 
Gaussian with common mean for AE. The CF and WC PDFs consist of sums of Gaussians 
and a Chebyshev function for AE, and Gaussians and either a Novosibirsk function [nl or a 

for M bc . The fractions of WC and CF-down events are taken from 



Crystal Ball function 



MC. The fractions of CF-up events are floated as they are difficult to simulate accurately 
(i.e., many B° s partial widths are unmeasured). As the CF-down fractions are fixed, the 
separate DfD~, D^DJ, and D* s + D*~ samples are fitted simultaneously. 

The projections of the fit are shown in Fig. and the fitted signal yields are listed in 
Table [H The branching fraction for channel % is calculated as B i = YJ {s' l MC -N B ^ ■ -2), 
where Y i is the fitted CR yield, and e\. IC is the MC signal efficiency with intermediate 
branching fractions [ijj included. The efficiencies e\. IC include small correction factors to 
account for differences between MC and data for kaon identification. Inserting all values 
gives the branc hing fractions listed in Table [H These results are in agreement with theoretical 



predictions [l9|, |20j. The statistical significance is calculated as J — 2 ln(£o/£ max ), where Lq 



and jC max are the values of the likelihood function when the signal yield Y i is fixed to zero 
and when it is floated, respectively. We include systematic uncertainties (discussed below) 
in the significance by smearing the likelihood function by a Gaussian having a width equal 
to the total systematic error related to the signal yield. 

The systematic errors are listed in Table [TT1 The error due to PDF shapes is evaluated by 
varying shape parameters by ±lcr. The errors for the fixed WC and CF-down fractions are 
evaluated by repeating the fit with each fixed fraction varied by ±20%. Those fractions that 
are correlated (e.g., WC for D*D+ and D* + D*~ , which is due to reconstructing extraneous 
photons) are varied together in the ratio predicted from MC simulation. The systematic 
errors due to qq suppression and the best candidate selection are evaluated using control 
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samples of B® — > D~n + and B° — > D^ + D~, respectively. These errors are taken as the 
change in the branching fractions when the criteria are applied. The uncertainties due to 
it ± /K ± identification and tracking efficiency are obtained from D* + — > D°ix + — > K~ir + 7T + 
decays; these are ~ 1% and 0.35% per track, respectively. Significant uncertainties arise 
from the T(5S) — > B*B* and branching fractions, which are external factors. We take 
the D* + D*~ polarization f L for this measurement to be the well-measured value from the 
analogous decay P^— >■ D* + D*~: 0.52± 0.05 jl^]. The systematic error is taken as the change 
in B when f L is varied over a wide range: from 2a higher than 0.52 down to the (low) central 
value we measure below. 

TABLE I: Signal yield (Y), efficiency including intermediate branching fractions (e), branching 
fraction (B), and signal significance (S) including systematic uncertainty. The first errors listed 
are statistical; the others are systematic. The last error for the sum is due to external factors 
(T(5S)^B*B* and D+ branching fractions). 



Mode 


Y 

(events) 


£ 

(xl0~ 4 ) 


B 

(%) 


5 


DfD- 


oq 1+6.0 
°-- I --5.4 


4.72 


0.58l£;o9 ±0.13 


11.5 


D*±DJ 


44.5iIJ 


2.08 


1.8 ±0.2 ±0.4 


10.1 


D* S D* S 


24.41" 


1.01 


2.0 ±0.3 ±0.5 


7.8 


Sum 


102.0±1;1 




4.3 ± 0.4 ±g;| ±0.9 





In the limits m, b c \ — > oo with {m b — 2m c ) — > and ^(number of colors) — > oo, the D^Df 
and DT^_D*~ modes are CP-even and (along with DfD~) saturate the width difference 
Arf p [2| . Assuming negligible CP violation (0 S ^O), the branching fraction is related to 
Ar s via Ar s /T s = 2i3 /( 1 — B). Inserting the total B from Table |T] gives 

AI\ 

— - = 0.090 ±0.009 ±0.023, (1) 

s 

where the first error is statistical and the second is systematic. This result has similar 
precision as that of other recent measurements 0, [21]]. The central value is consistent with, 
but lower than, the theoretical prediction (4]; the difference may be due to the unknown CP- 
odd component in B° s — >D* + D*~, and contributions from three-body final states. With more 
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FIG. 1: M hc and AE projections of the fit result. The top row shows B® — > DfD s ; the middle 
row shows Bg—tDg^D^' and the bottom row shows B^D* S + D*-. The red dashed curves show 
CR+WC signal; the blue and purple solid curves show CF; the grey solid curves show background; 
and the black solid curves show the total. 

data these unknowns can be measured. The former is estimated to be only 6% for analogous 
B° -^D* + D*~ decays |22|], but the latter can be significant: Ref. j2fj calculates AT(B S -> 
£(*)£>(*) If (*))/r s = 0.064 ± 0.047. This calculation predicts Ar,/r s from £>(*)+£>«- alone 
to be 0.102 ± 0.030, which agrees well with our result. If the CP- violating phase <p s is 
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TABLE II: Systematic errors (%). Those listed in the top section affect the signal yield and thus 
the signal significance. 



Source 






D* S D 


s 


D*i 


-D*- 




+a 


—a 


+a 


—a 


+ CT 


—a 


Signal PDF shape 


2.7 


2.2 


2.2 


2.4 


5.1 


3.8 


Bckgrnd PDF shape 


1.5 


1.3 


1.3 


1.4 


2.9 


2.8 


WC + CF fraction 


0.5 


0.5 


4.7 


4.5 


11.0 


9.7 


qq suppression 


3.1 


0.0 


0.0 


2.7 


0.0 


2.1 


Best cand. selection 


5.5 


0.0 


1.5 


0.0 


1.5 


0.0 


ir^/K* identif. 


7.0 


7.0 


7.0 


7.0 


7.0 


7.0 


K$ reconstruction 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


7T° reconstruction 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


7 






3.8 


3.8 


7.6 


7.6 


Tracking 


2.2 


2.2 


2.2 


2.2 


2.2 


2.2 


Polarization 


0.0 


0.0 


0.8 


2.4 


0.4 


0.2 


MC statistics for e 


0.2 


0.2 


0.4 


0.4 


0.5 


0.5 


br. fractions 


8.6 


8.6 


8.6 


8.6 


8.7 


8.7 


f ~* 






18.3 
2.0 








Total 


22.7 


21.8 


22.7 


22.9 


26.2 


25.5 



nonzero, Ar s /r s = 1/ cos0 s — J (1/ cos0 s ) 2 — 413(1 — B)/(l — B), where B includes all 
CP-even decay modes. Figure [2] plots Ar s as a function of <j) s for our measured B value. 

To measure f L , we select events using the same criteria as before but, to minimize 
B® — > D* DJ cross-feed, we use a narrower range of M bc and AE (2.5a in resolu- 
tion). For these events we perform an unbinned ML fit to the heiicity angles 9 l and 
9 2 , which are the angles between the daughter 7 momentum and the opposite of the B s 
momentum in the D* + and D*~ rest frames, respectively. The angular distribution is 
(\A + \ 2 + |A„| 2 ) (cos 2 6>i + 1) (cos 2 6» 2 + 1) + \A Q \ 2 4sm 2 9 1 sm 2 9 2 , where A + , A_, and A are 
the three polarization amplitudes in the heiicity basis. The fraction f L = |ylo| 2 /(| A)| 2 + 
\A + \ 2 + |t4_| 2 ). To account for resolution and efficiency variation, the signal PDFs are taken 
from MC. The background PDF is taken from an M bc sideband; the level (1.8 ± 0.7 events) 
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FIG. 2: The width difference Ar s as a function of <f> s . The dotted contours represent ±1<7 
bands from our measurement. We also plot the SM theory value [J], and the theory value for 
B s ^Di* )+ Di* ] ~ only 



20T | . The former is offset from the tiny SM 4> s value to make the full band 



visible. 



is estimated from a Df mass sideband and fixed in the fit. We obtain 



fl 



n 1 9 +0.30 +0.04 
u - iz -0.28 -0.05 ' 



(2) 



where the systematic errors arise from: signal PDF shapes (+0.016,-0.013), background 
PDF shape (±0.018), fixed WC fraction (+0.009, -0.011), fixed background level (±0.028), 
qq suppression (+0, —0.027), and MC efficiency due to statistics (±0.001). The helicity angle 
distributions and fit projections are shown in Fig. EJ 

In summary, we have measured the branching fractions for B® — yDy) + Dy' using e + e~ 
data taken at the T(55*) resonance. Under some theoretical assumptions and neglecting CP 
violation, the total branching fraction gives a constraint on Ar s /r s . We have also made the 
first measurement of the B® — > D* + D*~ longitudinal polarization fraction. 
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FIG. 3: Helicity angle distributions and projections of the fit result. The green dashed (red 
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